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ABSTRACT The relation of combinants to various statistics 
characterizing the fluctuation pattern of multihadron final states 
is discussed. Scaling laws are derived for count probabilities 
and combinants in the presence of homogeneous and clustered 
monofractal density fiuctuations. It is argued that both types 
of scaling rules are well suited to signal Quark-Gluon Plasma 
formation in a second-order QCD phase transition. 
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1. Introduction 

There is a renewed interest in an almost forgotten statistic well suited to 
characterize the fluctuation pattern and the nature of correlations of multihadron 
final states. This statistic is the set of combinants introduced 15 years ago by 
Gyulassy and Kauffmann [1-2]. After their introduction the combinants received 
little attention, only three papers have been published concerning these quantities 
[3-5]. Probably this was caused by the fact that the existence of combinants 
requires nonzero probability for detecting no particles at all (and precise data for 
this probability) which was the exception rather than the rule 15 years ago in 
the analysis of multiplicity distributions. It is worth mentioning that the Buras- 
expansion of multiplicity distributions invented years earlier [6] also utilizes the 
combinants but the emphasis in this work was made on the factorial cumulant 
moments and the above important requirement was unnoticed. This expansion was 
applied in [7-8] . Three further publications investigating the combinants have been 
appeared quite recently [9-11]. The work of Szapudi and Szalay [10] is particularly 
recommended to the interested reader. Although not explicitly mentioned the 
combinants are present also in the work of Lupia, Giovannini and Ugoccioni [12]. 

This letter concentrates on the study of various scaling laws of count prob- 
abilities and combinants that appear in the presence of monofractal density fluc- 
tuations. The discussion of the relationship between combinants and the more 
conventional measures of fluctuations and correlations is also the goal of the pa- 
per. Section 2 summarizes the basic properties of the combinants and a scaling 
law is provided which follows from the validity of the Linked Pair Approximation 
(LPA). A similar scaling law can be derived for count probabilities if the observed 
events have monofractal structure. This result will be presented in section 3. In 
section 4 the two types of scaling rules are reinvestigated for a Poisson superposi- 
tion of clusters. Our conclusions are in section 5. 

2. Combinants, LPA and scaling 

The generating function Q{z) of the multiplicity distribution P„ is conve- 
niently defined by 

oo 

g{z) = Y,PnZ^. (2.1) 

n=0 

In the analysis of various count probability distributions occuring in nature such as 
the distribution of galaxy-, hadron- and photocounts the most successfully applied 
distribution functions are infinitely divisible (e.g. Poisson, Gaussian, Negative 
Binomial, Lognormal). For discrete distributions this family is known also as the 
compound Poisson [13]. A discrete distribution is said to be infinitely divisible if 
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its generating function has the property that for all /c > integer y^Q{z) is again 
the generating function of a certain distribution [13]. Q{z) satisfies this property 
if and only if Q{1) = 1 and 

oo 

iiig{z)^iiig{o) + J2Cqz'i (2-2) 



where Cq > and 



oo 

E 

q=l 



Cg = -ln^(O) < OO. (2.3) 



Hence the probability of detecting no particles, the so-called void probability, is 
Pq = > for infinitely divisible distributions. The quantities Cq in eqs. 

(2.2-3) are the combinants first studied in detail by Gyulassy and Kauffmann 
[1-2] . We can get the Cq by successively differentiating the logarithmic generating 
function: 

and Q{z) takes the form 

g{z)=e^p(f^Cq{z'^-l)] (2.5) 

for infinitely divisible distributions. It is worth to introduce the count probability 
ratios Vn = {Pn/Po)- These can be expressed in terms of combinants by the 
recursion [10] 

1 

Vn^-Y^qCqVn-q. (2.6) 
q=l 

The Vn involve only a finite number of combinants but Pq (hence each P„) requires 
all the Cq as is seen from eq. (2.3). In terms of probabilities the combinants are 
found to be 

1 "'^ 

Cq = Vq--y^nCnPq-n- (2-7) 

Essentially this expression was obtained also in ref. [12] in a different context. From 
eq. (2.7) two advantageous features of the combinants are immediately seen. First, 
they require the knowledge of only a finite number of probabilities Pn- In g-th 
order the P<q are required. Second, we need not know the probabilities themselves. 
The combinants follow directly from the unnormalized topological cross sections 
since they involve only ratios of probabilities. It is also seen from eq. (2.7) that 
in general the combinants can take negative values as well for q > 2. In this case 
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a necessary condition of the infinite divisibility of the multiphcity distribution 
is not satisfied. 

In the analysis of correlations in restricted domains of phase-space a basic 
set of quantities are the integrated cumulant correlation functions providing the 
factorial cumulant moments, /q, of the underlying multiplicity distribution [14-15]. 
The factorial cumulants are defined by the Taylor expansion of the logarithmic 
generating function: 

l^^(^) = E^^/«- (2-8) 

q=l 

One of the main advantages of the combinants is their close relationship to the 
factorial cumulant moments. This manifests itself e.g. in the very simple expres- 
sions connecting the two quantities. The factorial cumulants can be expressed in 
terms of combinants according to 

(2.9) 

n=q 

and the combinants read as 

q=n ^ ^ 

in terms of factorial cumulants [2,6,10]. They also share some important common 
features. For the Poisson distribution we have Ci = n and C>2 = 0. Hence nonva- 
nishing higher-order combinants measure the degree of deviation from Poissonian 
behaviour. The combinants exhibit the additivity property of the cumulant mo- 
ments. When a random variable is composed of statistically independent random 
variables the corresponding combinants are additive in those of the independent 
components [2]. 

In the last few years much interest has been devoted to the analysis of correla- 
tions in the framework of the Linked Pair Approximation [16-17]. In the LPA the 
normalized factorial cumulant moments, Kq — fq/n'^, obey the recurrence relation 

Kq = AgRf^ (2.11) 

where the coefficients Aq should ideally be independent of reaction type, energy, 
binsize and phase-space dimension. In hadron-hadron collisions at ^/s = 200 - 900 
GeV the validity of LPA was confirmed with constant ^3^4 over the investigated 
range of energy and binsize [18]. The large uncertainties of A5 permit a clear 
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conclusion. At lower energies, = 22 GeV, the constancy of the Ag over binsizes 
was found to be violated by the NA22 Collaboration [19]. This may be caused 
by the absence of translation invariance. In nucleus-nucleus collisions the Ag are 
essentially zero indicating that only two-particle correlations are present [20] . 

Beyond the direct test of eq. (2.11) using the factorial cumulants there is 
another possibility to check the LPA scheme. It is provided by a scaling feature 
of the combinants [11]. As is seen from eq. (2.10) if the Kg obey the LPA relation 
the ratios Cn/n chosen at fixed n should scale to a function of nK2 only: 



Increasing n the contribution of an increasing number of low-order factorial cu- 
mulant moments is excluded from the Xn- Eq. (2.12) provides the generalization 
of the scaling law obtained for the void probability in the framework of the LPA 
[21,12]. As we shall see scaling combinants may appear also in the presence of 
monofractal density fluctuations. 

3. Monofractals and scaling count probabilities 

In recent years the observations revealed that in all types of reactions large 
nonstatistical fluctuations exist in the local particle density [22]. Varying the 
resolving power, i.e. the size of the phase-space domain in which the particles 
are counted the density irregularities resembled the scale-invariant, intermittent 
fluctuations that appear in a turbulent fluid [23-24]. This manifests itself in the 
power-law dependence of the normalized factorial moments on the resolution char- 
acterized by the anomalous dimensions dg or intermittency exponents {q — l)dq. 
The factorial moments, ^g, are defined by the Taylor expansion of Q{z): 



Measuring factorial moments in a certain phase-space domain is equivalent to 
measuring the underlying multiplicity distribution. Let us consider the relation 
between the $,g and the Pn- The factorial moments are expressed in terms of 
probabilities according to 



where the summation produces the binomial moments. The count probabilities 
take the form 




(2.12) 




(3.1) 




(3.2) 




(3.3) 
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when expressed in terms of factorial moments. Comparing eqs. (3.2-3) to eqs. 
(2.9-10) we see that the relationship between probabilities and factorial moments, 
generated by G{z) and Q{l+z), is the same as the relationship between combinants 
and factorial cumulants, generated by lnQ{z) and ln^(l + z). This important fea- 
ture which is valid for ordinary moments and cumulants too was first emphasized 
by Kauffmann and Gyulassy [2] . 

The above correspondence enables us to apply the same scaling argument for 
count probabilities which was utilized in the previous section for combinants. Let 
us assume that instead of the normalized factorial cumulants, Kg, the normalized 
factorial moments Fq = ^q/n'^ obey the LPA-type recurrence relation with constant 
coefficients Aq-. 

Fq = AqFr\ (3.4) 

Eq. (3.4) is a special case of the famous Ochs-Wosiek empirical relation [25] that 
can be formulated as 

Fq=AqF^''~^^'^''''^\ (3.5) 

From currently available data for the factorial moments, Aq ^ 1 and the inter- 
mittency exponent ratios {q — l)dq/d2 are largely independent of reaction type, 
energy and phase-space dimension [22]. If F2 changes according to a power-law 
with varying resolution eq. (3.4) corresponds to monofractal density fluctuations 
characterized by dq/d2 = 1, i.e. by the unique anomalous dimension d2 given 
by the intermittency exponent of F2. Plugging into eq. (3.3) with ^q expressed 
according to eq. (3.4) and evaluating the ratios Pn/n one arrives at 

^ = E f ^) ^-^^Aq imr-' = Vn{nF2). (3.6) 

q—n ^ ^ 

That is to say, if eq. (3.4) holds for the factorial moments the Pn/n chosen at fixed 
n not depend on reaction type, energy, binsize and dimension arbitrarily, but only 
through the momentum combination nF2. Picking up the Pn from the tail of the 
multiplicity distributions one subtracts the contribution of the low-order factorial 
moments from the scaling functions. In this manner we are able to go further in 
testing the monofractality of density fiuctuations than the highest order factorial 
moment that can be extracted from observational data. Since monofractal patterns 
in multihadron final states are expected to form during a second-order QCD phase 
transition the collapse of the Pn/n to the scaling curves rjn{nF2) as dq/d2 1 
may serve as a signature of Quark-Gluon Plasma formation. This point will be 
discussed in the next section. We note that the absence of translation invariance 
can lead to the violation of the binsize-independence of Aq and thus the scaling 
behaviour of P„- 
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4. Scaling laws in the Poisson cluster model 

Let us assume that the observed events are composed of identical groups of 
particles (clusters, clans, Quark-Gluon Plasma droplets or whatever) distributed 
according to a Poisson process. The generating function of the total event multi- 
plicity distribution P„ becomes [26-27] 

g{z) = exp {C{n{z) - 1)) (4.1) 

which is the convolution of the Poissonian generating function of the distribution 
of clusters having mean C and the generating function ^{{z) of the distribution of 
particles within a single cluster. Eq. (4.1) is another way of writing the generating 
function of infinitely divisible distributions, eq. (2.5), with C = — In^(O) and 

In eq. (4.2) the g-particle count probability in a cluster, pq, is found to be [10-11] 

Pq - (4.3) 

Since ?i(0) = each cluster must contain at least one particle, i.e. po = Co = 0. 
The two basic parameters of the model, the average cluster multiplicity, C, and 
the average multiplicity in a cluster, q, are given by [28,12] 

C = - InPo and q = n/C. (4.4) 

In the Poisson cluster picture the normalized factorial moments of a single cluster, 
Tq, can be expressed in terms of C and the normalized factorial cumulants Kq of 
the total events as [29-30] 

J^q = C'-^Kq. (4.5) 

Writing Kq according to the LPA relation, eq. (2.11), yields [28,11] 

Tq = AqJ^r\ (4.6) 

Thus we see that a Poisson superposition of clusters which obey the Ochs-Wosiek 
relation, eq. (3.5), with intermittency exponent ratios {q — 1) leads to the validity 
of the LPA relation, eq. (2.11), for the total events. Comparing eqs. (2.3) and 
(4.3-5) one finds the following relationship between the scaling functions derived 
in the previous sections: 



(4.7) 
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The two sets of scaling functions appearing on the two different levels are equiva- 
lent. On the basis of cqs. (4.5-7) let us collect the correspondences between total 
event observables and quantities characterizing a single cluster: 

If the events are Poisson superpositions of self similar clusters we observe the 
power-law dependence of Kg instead of Fg (the Fg show a bending upward 
behaviour on log- log plot). Through the slopes of the Kg we actually see the 
slopes of the J-'g, that is, the fractal structure of a single cluster [29-30]. 

a) When the Poisson clusters underlying the events have monofractal structure 
and the J^g obey eq. (3.4) the Kg exhibit the LPA relation, eq. (2.11). Thus 
the validity of LPA may arise as the consequence of randomly superimposed 
monofractal clusters [28,11]. 

Hi) Accordingly, by the observation of Xg-scaling for the total events confirming 
the validity of LPA one actually sees r/g-scaling for single clusters confirming 
their monofractal structure. 

The search of monofractal patterns in multihadron final states plays a distinguished 
role in intermittency analyses. It has been conjectured that monofractal density 
fiuctuations could signal a phase transition from the Quark-Gluon Plasma [31-32]. 
At the critical point of a second-order QGP phase transition to the hadron gas in- 
termittent fiuctuations are expected to occur characterized by a unique anomalous 
dimension as in the Ising model [33]. Bialas and Hwa proposed the experimental 
confirmation of the validity of eq. (3.4) as a QGP signature, i.e. the g-independence 
of the anomalous dimensions dg in eq. (3.5) [31]. Another proposal, put forward by 
Peschanski, is based on the description of the monofractal fiuctuations at the criti- 
cal point in terms of clusters [32] . The unique anomalous dimension is interpreted 
in this case as the fractal dimension of the random set of intermittent clusters of 
the ordered phase inside the disordered one at the transition. For each proposal 
one has to assume that the monofractal density fiuctuations survive the further 
evolution of the system toward the hadronic final state. According to eq. (4.7) 
both sets of scaling functions could signal QGP formation. Events possessing ho- 
mogeneous monofractal fiuctuations exhibit ?7g-scaling whereas events composed 
of random monofractal clusters give rise (through the ?7g-scaling of the clusters) to 
Xg-scaling for the observed events. 

5. Conclusions 

Over the last few years the analysis of fiuctuations and correlations in re- 
stricted domains of phase-space revealed that the higher-order statistics charac- 
terizing the fiuctuation pattern are frequently related to two-particle statistics in 
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a very simple manner [22]. One example is the LPA relation, eq. (2.11), for the 
factorial cumulant moments. The validity of LPA was confirmed for pp collisions 
at CERN collider energies with constant coefficients Ag close in magnitude to the 
Negative Binomial values Aq — {q — 1)\ [18]. Another frequently encountered re- 
lationship between second- and higher-order statistics is the Ochs-Wosiek relation 
for the factorial moments. In addition to the simplicity of eq. (3.5) it has a kind 
of universality: rather independently of reaction type, energy, binsize and phase- 
space dimension, the coefficients Aq ^ 1 and the intermittency exponent ratios 
(g — l)dq/d2 are compatible with a Levy-stable law with Levy-index = 1.6 [22]. 

Eq. (3.4), the Ochs-Wosiek relation with intermittency exponent ratios {q — 1) 
characterizes monofractal density fluctuations. Such patterns are expected to 
form during a second-order phase transition from the Quark-Gluon Plasma to 
the hadron gas [31-32]. Recent studies of the phase structure of QCD suggest that 
for two massless quark flavours QCD undergoes a second-order transition and this 
remains valid for three flavours provided that the strange quark mass is sufflciently 
large [34]. Hence it is of interest to find clear signatures of monofractal density 
fluctuations appearing in multihadron final states. 

In this paper we have derived a scaling law for the count probabilities P„ in 
the presence of monofractal density fluctuations. It is expressed by eq. (3.6). Ac- 
cording to this scaling law, if one plots the ratios Pn/n chosen at flxed n against the 
momentum combination nF2 the validity of eq. (3.4) with constant coefficients Aq 
results in a universal curve, rjn{nF2), instead of many different behaviours corre- 
sponding to different reaction types, energies, binsizes and phase-space dimensions. 
Increasing n the low-order factorial moments can be excluded in a systematic man- 
ner from testing eq. (3.4). We have also found that a similar scaling law is satisfied 
for monofractal density fluctuations occuring in randomly distributed clusters. In 
this case the combinants show up a scaling behaviour: at flxed n the ratios Cn/n 
fall onto a universal curve when plotted against the momentum combination nK2- 
Through the scaling curves one observes the validity of eqs. (3.4) and (3.6) again 
but now for a single cluster. Both types of scaling rules provide easily measurable, 
clear signatures of homogeneous and clustered monofractal density fluctuations 
which could be the remnants of a second-order QCD phase transition. It will 
be interesting to see how these scaling predictions are satisfled in next heavy ion 
experiments. 
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